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5-H ■ We propose a full-rate iterated space-time code construction, to design codes of Q-rank 2n from cyclic algebra 

based codes of Q-rank n. We give a condition for determining whether the resulting codes satisfy the full diversity 
property, and study their maximum likelihood decoding complexity with respect to sphere decoding. Particular 
emphasis is given to the asymmetric MIDO (multiple input double output) codes. In the process, we derive an 
interesting way of obtaining division algebras, and study their center and maximal subfield. 

Index Terms 

Cyclic Algebras, Division Algebras, Fast Decodability, Full Diversity, Quaternion Algebras, Space-Time 
O ' Coding. 



I. Introduction 

We consider the problem arising in wireless communication of coding for a multiple antenna system, where the 
transmitter has twice more antennas than the receiver. Codes for that setting fall into the category of asymmetric 

in 

space-time codes, and are usually formed of 2n x 2n complex matrices. The case when the receiver has exactly 
two antennas is called a MIDO channel (which stands for multiple input double output). Particular attention has 
recently been paid to the MIDO case, due to its potential application to digital TV broadcast, where the end user 
| carries a portable TV device. 

Many good space-time codes are available in the literature, in that they satisfy the full diversity property requiring 
that dct(X — X') ^ for any two matrices X ^ X' in the codebook - ensuring good behavior at high signal-to- 
noise ratio (SNR), and the non-vanishing determinant (NVD) property, which states that minx^x' | det(X — X')\ 
is bounded away from zero, even for an infinite space-time code, thus guaranteeing that the code performs well 
irrespectively of the alphabet size. 

However, most of the codes enjoying the above properties in fact possess an underlying lattice structure, and might 
suffer from high maximum likelihood (ML) decoding complexity. This triggered a new line of research focusing on 
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Communication, Control, and Computing, 2011, Illinois, USA, and in the proceedings of the International Conference on Signal Processing and 
Communications (SPCOM 2012), Bangalore, India. 
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designing codes with both good performance and reduced ML decoding complexity via sphere decoding [ l j, called 
fast decodable codes J2). Analysis of conditions under which a code satisfies fast decodability can be found in 
e -g- 0-GI, where the notions of fast decodable, group decodable and conditionally group decodable are detailed. 
Several 4-dimensional fast decodable MIDO codes have subsequently been studied, attempting to improve on the 
MIDO code in \2\, which combines a quasi-orthogonal code with a twisted unitary transformation of another quasi- 
orthogonal code, and gets a lower decoding complexity, by sacrificing the full diversity property, and consequently 
the stronger NVD property. In JS), cyclic algebra based codewords are arranged in an Alamouti block code J7] 
fashion to gain fast decodability for asymmetric space-time codes. Recent fast decodable code constructions include 
a fully diverse MIDO code conjectured to have the NVD property (8|, and MIDO codes with the NVD property 
based on crossed product algebras |0, iTTOl . Constructions for arbitrary number of transmit antennas are also 
available in iTTTl . 



A. System Model and Fast Decodability 

We consider transmission over a coherent Rayleigh fading channel with 2n transmit antennas, n receive antennas 
and perfect channel state information at the receiver (CSIR): 

Y nX 2n = ^11X211^211X211 + VnX2ni (1) 

where H is the channel matrix and V is the noise at the receiver, both with i.i.d zero mean complex Gaussian 
components. The 2n x 2?? space-time code X can transmit up to 2n 2 complex (say QAM) information symbols, or 
equivalently An 2 real (say PAM) information symbols. We will focus on the full-rate case where indeed An 2 real 
symbols are sent. When n = 2, X is sometimes referred to as MIDO space-time code, where MIDO stands for 
multiple input double output. 

Maximum-likelihood (ML) decoding consists of finding the codeword X that achieves the minimum of the 
squared Frobenius norm 

d(X) = ||l^ix2n — H n x2nX\\ 2 F . (2) 

By describing the space-time code X in terms of basis matrices i = I,..., An 2 , and a PAM vector g = 

(9i,---,9in 2 ) T as 

4ra 2 

X = Y,9iBi, (3) 

i=l 

we can rewrite (f2) in Euclidean norm as 

d(X) = ||y-5g||i, (4) 

where y <S K 4 " is the channel output which has been vectorized, with real and imaginary parts separated, and 

B = (bx.ba, • ■ • ,b 4 „ 2 ) g M 4 „2 X4 „2(M), (5) 

also obtained by vectorizing and separating the real and imaginary parts of H nX 2 n Bi to obtain b;, i = 1, . . . , An 2 . 
This search can then be performed using a real sphere decoder (TJ. 
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We will focus on two aspects of space-time code design: diversity and fast decodability. 

Full diversity. It is well known lfl2l that the first important space-time code design criterion is full diversity, 
that is, we require det(X — X') ^ 0, for any X ^ X' in the codebook. When X ± X' is again a codeword for all 
X, X', the full diversity criterion simplifies to 

det(X) ^ for all X ^ 0. (6) 

Fast decodability. From @, a QR decomposition of B, B = QR, with Q unitary, reduces to computing 

d{X) = | |y - QRg\ || = ||Q*y - Rg\ || (7) 

where R is an upper right triangular matrix and ()* denotes the Hermitian transpose. The number and position of 
nonzero elements in the upper right part of R will determine the complexity of the sphere decoding process. 

Definition 1: If S is the real alphabet in use, and k is the number of independent real information symbols from 
S within one code matrix, then the ML decoding complexity [|2[ is the minimum number of values of d(X) in (0 
that should be computed while performing ML decoding. 

Note that n will sometimes be referred to as the rank or dimension of the code, which corresponds to the number 
of basis matrices, k = 4n 2 in our case. The worst case is given when the matrix R is a full upper right triangular 
matrix, yielding the complexity of the exhaustive-search ML decoder, that is here Od^l 4 ™ 2 ), with k = An 2 and 
| S | is the number of PAM symbols in use. If the structure of the code is such that the decoding complexity has an 
exponent of 15*1 smaller than An 2 , we say that the code is fast-decodable. 

We now recall three shapes of matrix R which have been shown O, H to result in reduced decoding complexity. 

Suppose a space-time code has matrix R of the form: 

A B 1 

R = 

R 2 

where A is a d x d diagonal matrix and R 2 is upper-triangular. Then decoding complexity of the code is reduced 
to 0(\S\ K ~ d+1 ). The matrix R, which depends on the channel matrix H, may be difficult to compute. Thankfully, 
the zero structure of R is actually related J3] to the zero structure of the matrix M defined by 

M kll = \\B k Bf+B l Bi\\F, (8) 

which captures information about orthogonality relations of the basis elements Bi, with B* denoting Hermitian 
transpose. It furthermore has the advantage that its zero structure is stable under premultiplication of Bi by a 
channel matrix H (in general, the same does not hold for R). It was shown in (3] Lemma 2] that 



A 






A 






R = 






B 2 


B 3 _ 







Ri 



where A is a d-dimensional diagonal matrix. It is thus sufficient to compute M to deduce fast-decodability. 

Definition 2: A space-time code of dimension k is called g-group decodable if there exists a partition of 
{!,...,«;} into g nonempty subsets Ti, . . . , T g , so that M k ,i = when k, I are in disjoint subsets I\, Tj. 
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In this case, as shown in |3), the matrix R has the form R = diag(Ri, . . . , R g ), where each Ri is a square upper 
triangular matrix. Hence, the symbols and gi can be decoded independently when their corresponding basis 
matrices Bk and B\ belong to disjoint subsets of the induced partition of the basis, yielding a decoding complexity 

of 0(|S , | maXl ^< s r -). 

Finally, a code might not be a g-group decodable code, but could become one assuming that a set of symbols 
are already decoded. This corresponds to another characterization of fast-decodability J4): 

Definition 3: A code is called conditionally g-group decodable if there exists a partition of {1, . . . , k} into g + 1 
disjoint subsets Ti, . . . , T g , T c such that 

\\B[B* m + B m Bf\\ F = vier.Vmerj,^ j. 

In this case, the sphere decoding complexity order reduces to | 1 S'| rC + lnaXl < i 5s Ti . 

By the above discussion, in order to demonstrate fast-decodability, respectively (conditional) g-group decodability, 
it suffices to find an ordering on the basis elements Bi, which results in the desired zero structure of M . 

B. Contribution and Organization 

In this paper, we consider the problem of designing full-rate space-time codes for 2n transmit and n receive 
asymmetric MIMO channels. Our main contributions are: 

1) We propose a general code construction that maps two rank-rt algebraic space-time codes coming from a 
cyclic algebra to a rank-2n new code and give a criterion for full diversity. 

2) We obtain a novel way of constructing division algebras, and describe their center and one maximal subfield. 

3) We provide a decoding complexity analysis that permits to derive conditions under which the iterated code 
inherits fast-decodable properties of the underlying code. 

4) We construct iterated MIDO versions of the well known Golden |[T3l and Silver code IU~4l , as well as MIMO 
codes for the much less studied case of 6 transmit and 3 receive antennas. 

This paper is organized as follows. The main construction is explained in Section[n] together with its main properties 
- full diversity, algebra structure and fast decodability. We focus on the case of quaternion algebras in Section 
IrTfl where we reformulate full diversity in terms of quadratic forms. MIDO codes constructions are provided in 
Section |IV| where we construct iterated versions of well known quaternion based codes, such as the Silver and 
Golden code. Finally in Section [V] cyclic algebras of degree 3 are investigated. Simulation results illustrate the 
performance of the codes. 

II. An Iterated Code Construction 

Let F be a finite Galois extension of the field Q of rational numbers. Let K/F be a cyclic Galois extension of 
degree n with Galois group generated by a. Let T> := [Kj F, <r, 7) = if © eK © • • • © e n_1 K be the cyclic central 
simple algebra over F defined by the relations 

e n = 7 e F, 
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le = ecr(0 VZ £ K. 

As a vector space, T> is of dimension n 2 over its center Z(D) = F. The algebra T> is said to be of degree n 
over its center F. We say that T> is division when every nonzero element of T> is invertible. 

Left multiplication of elements of T> induces the representation A in M nxn (K), which maps elements of T> into 
the ring M nxn {K) of n x n matrices with entries in K 03, given by 

I 
a(l) 



A : / n- 



cr"- 1 ^) 



,Z G K 



A : e H> 





1 



1 

Any automorphism r of K extends naturally to an automorphism of M nxn (K), also denoted by r 

t : (Aij) h-> (r(^-)) 
by applying r on the matrix coefficients. Additionally assume the following: 



(9) 



Then the automorphism r of M nxn (K) restricts to an automorphism of X(T>). Alternatively, the automorphism r 
of K extends to the automorphism r of V: 

t : xo + exi H + e n ~ 1 x n -i 

t(x ) + er(xi) + . . . + e' l_1 r(x„-i) 

and hence t extends to the representation \(D) via r : \(z) h-> A(r(z)). Putting the above together we have the 
following commutative diagram 

V : X(V) M nXn (K) 



V 



\{V) 



M nxn (K). 



Note that for matrices A, B e A(X>), t(A)t(B) = t(AB). Hence for A invertible, we have r(A _1 ) = r(A) -1 . 
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A. The General Code Construction 

Given an element 9 of T> and a Q-automorphism r of K, let 

ag ■ M nxn {K) x M nxn (K) -> M 2n x2n(K) 

be the map defined by 

X 9t(Y) 

ag :(X,Y)^ 

Y r(X) 

If in addition the assumptions in (O hold, then the restriction of ag to \(D) x \(D) can be seen as an embedding 
of V x T> into M2nx2n(K) by identifying elements x, y G T> with their representation \(x) = X,X(y) = Y E 
M nxn (K). 

Remark 4: In fact, this construction can be made more general in two ways: (1) one could start with a noncyclic 
central simple algebra (this has been proposed in |[T6l ). and (2) the map ag could be applied iteratively to obtain 
a family of codes for 2' transmit antennas, i > 0. ifTTl . 

Remark 5: Note that the case when r ^ a and both are Galois automorphisms of K of order 2, the iterated 
construction coming from a$(T>,T>) is in fact the same as the matrix representation of the crossed product 
(K/F, 7, 9, 1) |[l"8l induced by left multiplication, where the automorphisms a, r generate the biquadratic extension 
K/F and V = (K/K^ , a, 7). 

We start with a simple but important criterion to decide when the image A of ag has the structure of a finite 
dimensional Q-algebra. 

Lemma 6: Let T> = (K/F, a, 7) and let r G Gal(K/F) be an automorphism such that assumptions in (0 hold, 
i.e., rer = err and r(7) = 7. Additionally, suppose r 2 = 1 and fix G F such that r(#) = 0. Then the image 
.4 = ae(X>, 2?) forms a Q-algebra. 

Proof: Note that the assumptions rcr = err and r(-f) = 7 allow us to extend the automorphism t from K 
to I?. Next observe that the image A of ag is both additively and multiplicatively closed. The additive closure is 
immediate: 

u 9t(v) 



ag{x,y) +ag(u,v 



x 9r(y) 
y t(x) 
ag(x + u,y + v). 



v t(u) 



For the multiplicative closure, we verify that 

xu + 9r(y)v x9t(v) + 0t(j/)t(u) 

a 9 (x,y)ag(u,v) = 

yu + t(x)v y9r(v) + t(x)t(u) 
= ag(xu + 9r(y)v, yu + t(x)v) 

using the assumptions that t(9) = 9, 9 G Z(D) = F and r 2 = 1. Hence A is a ring. Since A is additionally a 
vector space over ag(Q, 0), which is central in .A, identifying Q with a#(Q, 0), gives the ring A the structure of 
a Q-algebra of (finite) dimension [A : Q] = 2 [2? : F] [F : Q] = 2n 2 [F : Q], ■ 
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We now give a condition on 9 which guarantees that the image A of ag satisfies the full diversity property (|6), 
irrespectively of whether A is an algebra (more precisely, A might not be multiplicatively closed). 

Lemma 7: Let K/F be a cyclic extension of number fields and T> = (K/F,a,j) be a cyclic division algebra 
and fix 9 £ T>. Let r be a Galois Q-automorphism of K satisfying assumptions in ([9j, that is, t(j) = 7, rer = err. 
Every matrix in the image A of ag has determinant in F and is invertible if and only if 9 7^ zt(z) for all z £ V. 
If moreover t 2 = 1 and 9 £ F, then for any a$(x,y) £ *4, the determinant det(ag(x, y)) is contained in the 
subfield F^ of F fixed by r. 

Proof: Suppose that 9 7^ zr(z) for all z £ V. We will show that A is fully diverse, i.e., every nonzero element 
of A has nonzero determinant. 

x 9r(y) 

Consider a nonzero element ag(x, y) = , where the entries x, y are n x n matrices from the division 

y t(x) 

algebra T>. It is enough to demonstrate that the determinant of ag(x,y) is nonzero, since A is additively closed. 

If x = (resp. y = 0), the matrix is clearly invertible. Hence for x,y £ A(X>), we assume that x and y are both 

A B 

nonzero, and hence, as elements of a division algebra, both invertible. The determinant of a block matrix 

when A is invertible, is given by det(A) det(D — CA~ l B). Then 
we have 

x 9r(y) 



C D 



det 



= det(x) det(r(a;) — yx 9r(y)). 



(ID 



y t(x) 

Note that t(-j) = 7, together with the fact that r and a commute, imply that t(x) € A(2?). The assumption 9 E F 
is needed so that 9r(y) £ A(2?). It then follows that t(x) — yx~ 1 9r(y) is an element of A(2?), and hence invertible 
when nonzero: we have dct(r(x) — yx~ 1 9r(y)) 7^ t(x) — yx~ 1 9r(y) 7^ 0. It thus suffices to demonstrate 

t(x) — yx~ 1 9r(y) ^ 0. Since y is invertible (it is nonzero), and noting that for x £ A(2?), t(x^ 1 ) = t(x)~ 1 , the 
latter inequality is equivalent to xy^ 1 r(xy^ 1 ) 7^ 9. Letting z = xy^ 1 , by our assumption on 9, we have 9 7^ zr(z). 
Hence det(ag(x, y)) 7^ and is well known to belong to the center F of T>, since both x and t(x) — yx~ 1 9r(y) 
are elements of T>. 

Now assume additionally that t 2 = 1 and 9 £ F. We show that det(ag(x, y)) £ F ^ . Since by above remarks 
we already know that the determinant is in F, it suffices to show that it is fixed by r. Now when D is invertible, 



we can alternatively express the determinant of the block matrix 
the assumption that t(A) is invertible with thus get 

9r{y) 



A B 
C D 



as det(£>)det(A - BD~ l C). Using 



dot 



= det(r(.T)) det(x — 9r(y)T(x) 1 y) 



(12) 



V r(x) 

Equating formulas (fTTT > and (fT2l makes it clear that det(ag(x, y)) is fixed by r. Hence we have det(ag(x,y)) £ 
Conversely, suppose 9 = zt{z) for some z £ X(T>). We verify that dct(ag(z, I n )) = 0, where /„ is the 
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n-dimensional identity matrix. Indeed 

det(ag(z,I n )) = dct(z) det(r(z) — I n z~ 1 9) 

= det(z) det(r(z) — l n z~ x zt(z)) 
= det(» det(O) = 0. 

Hence A is not fully diverse. ■ 
Corollary 1: Let T> = (K/F, er, 7) be a cyclic division algebra. Consider a Galois Q-automorphism r of A', 
such that assumptions (O hold, i.e., r(~f) = 7, rcr = err. Assume additionally that t 2 = 1 and fix G F such that 
t(6>) = 9. Then the image A of is an algebra, which is division if and only if 9 ^ zt(z) for all z G T>. 

Proof: The assumptions ©onr allow to extend t to an automorphism of T>, while the additional assumptions 
that t 2 = 1,9 G F and t(9) = 9 give the image A the structure of a finite dimensional Q-algebra by Lemma [6] 
A finite dimensional algebra is division (i.e. every nonzero element has an inverse) if and only if it contains no 
zero divisors. If 9 ^ zt(z) for all z G T>, then by previous Lemma every element of A is an invertible matrix, and 
hence cannot be a zero divisor. Hence A is division (note that we avoided having to explicitly demonstrate that 
the inverse of every element is indeed in .4). On the other hand if 9 = zt(z) for some z G V, then by previous 
Lemma A contains a non-invertible matrix, and therefore cannot be division. ■ 

B. More Algebraic Properties and Fast Decodability 

We saw in Lemma [6] that assuming that r 2 = 1 and 9 G F< r >, the image A is an algebra. It is natural to wonder 
what is the center Z(A) of A. 

Lemma 8: Let r be a Galois automorphism of K of order 2. If r G (cr), then the ring ag(F, e n ^ 2 F) is a maximal 
subring which commutes with all elements of A. Otherwise, ag(F T ,0) is such a ring. 

Proof: The Galois automorphism r may or may not be an element of the Galois group (a). We treat the two 
cases separately. 

Case 1. Let us first suppose the former case when r G (cr). Since r is of order 2, it follows that r = a n l 2 . 

Suppose that ag(u,v) commutes with all the elements of A, that is, a$(u, v)ae(x, y) = ag(x,y)ae(u,v) for 
any x, y G T>. Then, in particular, it must be true for y = 0. This implies that ux = xu and t(x)v = vx must hold 
for all x G T>. The first condition implies that u G Z(V) = F. The second implies that a n ' 2 (x)v = vx. 

We show that this condition forces v to be an element in e n / 2 F. Express v = vq + ev\ + • • • + e n ~ 1 v n ~ 1 , where 
vi, . . . , v n _i G K. Now for all x G K, we have 

a n ' 2 {x)v = a n ' 2 {x){v + evi + ■ ■ ■ + e™" 1 ^!) = 

v (T n/2 {x) + ev x a n ' 2+1 {x) + ■■■ + e r ^ 1 w„_ia"/ 2 +"- 1 (x). 

Equating with vx = vqx + ev\x + ■ ■ ■ + e n ~ 1 v n _ix in each component gives e l ViO- n l 2+l (x) = e l ViX for all x G K. 
This forces for all vi 7^ the condition a n / 2+l {x) = x for all x G K. The latter is true if and only if i = n/2, 
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A 




7£f 



Fig. 1. Tower of extensions showing the iterated construction: on the left, T> is of dimension n 2 over its center F. In the middle, when 
T £ (cr), A is of dimension 2 over T>, that is of dimension 2ra 2 over F. It is further of dimension n 2 over the algebra a(F, e n / 2 F), e n = 7, 
which is itself of dimension 2 over F. On the right, when r (cr), the algebra A is of dimension 4n 2 over _F T . 

forcing ^ for all i ^ n/2 to be zero. Hence we obtain v £ e n l 2 K. Now letting x = e and writing v = e n ^ 2 v n / 2 , 
the constraint a n l 2 (e)v = ue, gives e ■ e n ^ 2 v n / 2 = e"^ 2 v n / 2 e = e "^ 2ecr ( u n/2X implying that u„/ 2 = cr ( v n/2)> i- e -> 
w £ e n / 2 F, as desired. On the other hand, we check that ag(F, e n l 2 F) in fact commutes with all elements of A. 
This completes the first case. 

Case 2. Now suppose that r does not belong to (a). If u) commutes with all the elements of A, the conclusion 
that u € F holds similarly to Case 1. 

Now we additionally need to impose the constraint t(x)v = vx for all x <G T>. By writing v = ^™T e* 1 '^ this 
yields for any a; e K the condition cr^-r^a;))^ = ViX for all i. Since crV is not identity on A" for any i, we 
conclude that = for all i, showing that v = 0. Next we check for which values of u the element ag(u,0) 
commutes with all elements ag(0, y) where y £ T>. This produces the constraint r(u)y = yu for all y £ T>, that 
is, it £ F T . On the other hand, ag(F T , 0) indeed commutes with all the elements of A, completing the proof. ■ 

Remark 9: In Corollary [TJ we derived a criterion which forces A to be division. When this criterion is satisfied, 
then in case r £ (cr), we have F' := ag(F,e n / 2 F) is necessarily a field and it corresponds to the center of A, 
while K' := a e {K,e n / 2 F) is a maximal subfield of A. In fact, we can see that A = (K' /F', cr, 7) and is generated 
as a vector space over K' by powers of the element a§(e, 0). In case when r g" (cr), the center of „4 is ae(F T , 0) 
and ag(K, 0) is a maximal commutative subfield of A The latter follows from the observation that ag(K, 0) is 
clearly commutative and its maximality follows from the fact that it is of correct degree 2n over the center. 
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The embedding ag : T> x T> — > M2 n (K) always has an additively closed image A, which forms a vector space 
over F of dimension 2n 2 . In fact, we can generate A over F by elements of the form 



{ag(u, 0),a e (0, u) : u = ui, . . .u n i £ B F (D)}, 
where ,Bf(2?) is an F-basis of the cyclic algebra T> = (K/F,a,j): 



(13) 



V = K © eK < 



where 



K, e" = 7 £ F 



*AiF 



A" = i^F © z/ 2 F © ■ 
by fixing an F-basis {y\, . . . , v n } of K. We can thus choose 

B F (T>) = j>i, . . . , v ni eft, ev n , e"~Vi, 



If F is a quadratic field extension, we can obtain a Q-basis Bq(T>) = {Dj, j = l,...,2n 2 } of T> using a 
Q-basis {ft.ft} of F: 

{Fx, ... , D n 2,D n 2 +1 , . . . , D 2 „ 2 } = /3iBj,(X>) U /3 2 S F (2?) 

{Hj8i, . . . , e"" 1 ^!,!/!^, . . . , e n " V&} 

When F is a field of higher degree, we let /?i , /3 2 be any two Q-independent elements of F, to define a Q-independent 
set 

{D 1 , . .. , D n 2,D n 2 +1 , . . . , D 2n 2} = 



{v 1 f3 1 ,...,e n 1 v n Pi,viP2 



This finally yields a generating set 



B(A) = {B j , j = l,...,4n 2 } = 

{ae(u,0), a s (0, u) :u = Di,.. .D 2n 2] 

and encoding as described in (01 consists of mapping a vector (g%, . . , , g^) of PAM information symbols (that is 
from Z) into a space-time codeword in A using the above basis matrices: 

2n 2 

■• - jflwO ^ -X = ^(flr»ae(A,0) + g 2n 2 +i ag(Q, Di)). 

i=l 

We will use the following scaling technique. Let u,v,x,y £ V. When 8 is either real or totally imaginary, we 
express 9 = (9' with 8' > 0, for some ( £ {±l,±i}. Then the map 

u (8't(v) 



a^gi : (u,v) i ^ 



v t(u) 
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and the map a^ygy denned by 

u (VO't(v) 
f&v t{u) 

satisfy the property that dct(a^e'(u, v)) = dct^a^y^u, v)), for all u,v. Note alternatively that 



a <V¥ : (u,v) i y 



1 


a e {u,v) 


1 






o Ve> 




o i/VW 



and in particular, the image of o^yg? retains the full diversity (resp. NVD) property of the image of a^e' 
Furthermore, assuming that complex conjugation commutes with r on elements of V, we get 



so that letting a denote 5^ v /g7 we have 
a(x,y)a(u,v)* -- 



u C,\/Wt(v) 
f&v t{u) 

x CVWr(y) 
V&y t(x) 



u* Vd'v* 
C*V¥t(v*) t(u*) 



CVPt(v*) t{u*) 



xu* + 6'r{yv*) ^f(F(xv* + (r(yu*)) 
9' (yu* + £*t(xv*)) / yv*+r(xu*) 
Since ( is a fourth root of unity, we notice that xv* + (r(yu*) = (((*xv* + r(yu*)) and 



y(x, y)a(u, v)* = a(xu* + & 'r(yv*),yu* + £*t(ot*)). 



This allows us to compute 



a(x, 0)a(u, 0)* = a(xu* , 0) ■ 
i(x,0)a(0,v)* = &(0,C*t(xv*)) = V¥ 
a(0,y)&(u,0y = a(0,yu*) = Ve 1 
a(0,y)a(0,v)* = a(9'r(yv*),0) = 9' 



xu* 
t(xu*) 

xv* 
C,*t{xv*) 

(r(yu*) 
yu* 

r(yv*) 
yv* 



(14) 



(15) 



(16) 



(17) 



(18) 



from which we see that fast decodable properties of the code coming from T> will be inherited by the image 
a^^iV^V). More precisely: 

Lemma 10: Let T> be a division algebra whose corresponding code has basis {Z?i, . . . , D 2n 2}. If 

D 3 D* k + D k D* = 

for some j, k, then the basis {B\, . . . , B^ n i\ of the code coming from A = ae(T>, T>) given by 

{B t = a(D u 0), B 2n 2 +t = 5(0, A), i = 1, . . . , 2n 2 } (19) 
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satisfies 

B 3 B k + B k B j = B 2n 2 +j B 2n2 + k + B 2n 2 +k B 2n 2 +j = °- 

Proof: Indeed, let {-Di, . . . , D 2n 2 } be a basis for the code T> and (fT9b be one for the code from A. We need 
to compute BjB^ + BkBj for two groups of indices: j,k S {1, . . . , 2n 2 } and j, k e {2n 2 + 1, . . . , 4n 2 }. In the 
first case, we use (T5[ to see that 



B j B k + B kB* 





r(D k D*) 

DjDl + D k D 








t{D 3 DI 




from which it is clear that the orthogonal relations among Dj determine those among Bj. The other equality is 
handled similarly, using ( TT~8T >. ■ 
Further improvements on these reductions in complexity may be possible for particular choices of £, r and are left 
for exploration in future work. 



III. The case of Quaternion Algebras 

In this section we focus on the case when the underlying division algebra is of degree 2, and can equivalently 
be viewed as a generalized quaternion algebra. Additionally, we will only consider the case when r = a. 

Let a,b 6 F x . We recall from lfl"9l that a (generalized) quaternion algebra over F, denoted by (a, 6)f, is a 
4-dimensional F- vector space with F-basis {1, i, j, k}, whose ring structure is determined by the rules 



a, j = b and ij 



-Ji 



(20) 



Let F be a number field and let Q = (a, *y)p be a generalized quaternion division algebra. We can alternatively view 
<5 as a cyclic division algebra (K/F, a, 7), where a maximal subfield of Q is a quadratic extension A" = F(y / a) 
of F with Galois group generated by the automorphism a : -y/a n- — -y/a. The diagram of Fig. [2] illustrates this 
tower of extensions. 

The fact that Q is a division algebra implies that 7 is not a norm of K fl9l . We can write 

Q = K ® jK, i 2 =7GF, 



and view elements of Q as 2 x 2 matrices over A via matrix representation: 

C 7 cr(d) 

c + ja i-> 

c? cr(c) 

The map as defined in ( [Tol l with r = cr is given here by 

a; 6»er(y) 



a e : (x,y) i-> 
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Fig. 2. Tower of extensions showing the iterated construction. On the left, Q is of dimension 4 over its center F. On the right, A = ctg(Q, Q) 
for t = a and 9 G F, is of dimension 2 over Q, that is of dimension 8 over F. It is further of dimension 4 over the subring a(F,jF), j 2 = 7, 
which is itself of dimension 2 over F. 



By identifying x and y with their matrix representation, ag can be seen as an embedding of Q x Q into M^lf). 
From Lemma [6] we have that A = ag(Q, Q) is an f-algebra whenever 9 E F. 
Example 1: Consider Q = (—1, — 1)q and write 

Q = Q(i)®jQ(i), f = -1. 

Note that —1 is not a norm of Q(i) over Q, i.e., — 1 ^ (/ + <?*)(/ — <?*)> for f,g&Q. Therefore Q is a division 
algebra; its maximal subfield Q(i) is an extension of degree 2 over its center Q with Galois group Gal(Q(i)/Q) 
cyclic of order 2. In this case, the generator a of the Galois group is the complex conjugation: a : x i-> x*, 

x G Q(i). 

We view the elements of Q as 2 x 2 matrices over Q(i) via matrix representation: 

c -d* 



(21) 



d c* 

We illustrate the iterated construction arising from the quaternion algebra Q = {— 1, — 1)q in Figure |3] where 9 E Q 
is chosen so that .4 = ag(Q, Q) is division, by Corollary [TJ The center and maximal subfield of „4 correspond to 
<%e(Q>jQ)> QeiQii), jQ(i))> respectively, as follows from Lemma [8] and Remark |9] 

A. First Criteria for Full Diversity 

The following lemma gives necessary and sufficient conditions on 9, 7 in order to assure that the image A of ag 
is a division algebra, following Corollary Q] Let K x denote the multiplicative group of nonzero elements in K. 
Lemma H: Let K = F(y/a), r y,6 G F. We have the following equivalence: 

9 7^ z(t(z) for any z G Q = (a, 7)f 
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A 



Q = (-1,-1) F a g (Q(i),ji 

2 

K = a 8 {Q,ji 

2 

F = 




Fig. 3. The iterated constraction from Example [T] when A is division. 



1) 6 ± det 



u ja(v) 
v cr(u) 



for any u,v E K such that v G \faF, and 



2) 9 7^ 7 (mod K x2 ), where (mod K x2 ) denotes taking the quotient by the subgroup K x2 of squares in K > 
Proof: Assume conditions 1) and 2). For the sake of contradiction, suppose there exists an element z e Q 



such that 9 = zcr(z), and write it as z 
This gives 



v c{u) 



, where u, v 6 K. 





9 




9 



ua(u) + 7<r(i') 2 U71> + ~fcr(v)u 
va(u) + a{u)a{v) v~/v + a(u)u 
Suppose « / 0. Comparing entries (1,2) of the above matrices (keeping in mind that u,v € K, and hence 
commute) gives v + a(v) = 0, which implies that v = vi^/a, for some t>i € F, since If is a quadratic extension 
of F . Comparing entries (1,1) gives 

ua(u) + jaiv 2 ) = 9. 
Equivalently, recalling the previous condition that v = V\^fa, we have 



9 = det 



u 7(7(11) 
v c(m) 



The first assumption on 9 gives us a contradiction. 

Now we consider the case u = 0. Comparing entries (1,1) gives us the equality 

7<t('u) 2 = 7<r(?; 2 ) = 9. 

By the second assumption that 9 7^ 7 (mod K x2 ), we obtain a contradiction. 



March 12, 2013 



DRAFT 



15 



For the converse, first suppose that 6 — ■ *yv for some element v £ K . We show that — where 

7ct(u 2 ) 



7<r(?;) 
u 



We have za{z) 



7« 2 



Note that since 9, 7 € F, the assumption that 



7?; 



implies that u 2 <E F. Hence v 2 = <r(v 2 ) and we have therefore just shown that 9 = za(z) 

u ja(v) 



Now suppose 9 = dot 



v c(u) 



where v G \faF. We already demonstrated how 9 = za(z), where z 



v cr(u) 

Note that since our construction allows the ground field F to be of arbitrary degree over (Q>, it can moreover be 
used for the multiblock construction in [20 1 . 

In the following corollary, we show cases where we can meet the conditions of the previous lemma. 
Corollary 2: Let F be a totally real number field, K = F(y/a), where a < 0. Let 7 < 0, and let 9 < satisfy 

9 7^ 7 (mod A' x2 ). Then nonzero elements of A have nonzero determinant. 

u 7(7(11) 

Proof: Note that by our choice of 7 < 0, a < 0, the determinant of is positive. Hence it cannot be 

v cr(u) 

equal to 9, which is by our choice negative. Together with the assumption that 9 ^ 7 (mod K x2 ), both conditions 
of Lemma QT| are satisfied and we conclude that nonzero elements of A have a nonzero determinant. ■ 
Example 2: In Example Q] F = Q, K = Q(i) and 7 = — 1. Thus any choice of 9 € F with 9 < not equal to 
— 1 (mod Q x2 ) will give a division algebra. 



B. Quadratic Forms Criteria for Full Diversity 

We rephrase the conditions of Lemma QT] in terms of quadratic forms. We can then use the theory of quadratic 
forms in order to meet these conditions. First we give some basic definitions. 

Definition 12: Let F be a field. A quadratic form on V over F is a map if : V — > F satisfying the following 
two properties: 

1) ip(av) = a 2 if(v) for all v e V and all a e F; 

2) the map b v : V x V — > F defined by btp(v, w) := <p(v + w) — f{v) — v( w ) f° r a U v,w 6 ^ is bilinear (and 
automatically symmetric). 

Definition 13: A quadratic form if on V is called anisotropic over F if it has no nontrivial solutions, i.e., 
if(v) = if and only if v = 0. 

We will only deal with quadratic forms which have the form ip(y) = Y^l=i a i v f f° r all v = (i>i, . . . , v n ) G F n 
and scalars 01, . . . , a„ € F, for convenience employing the short-hand notation ip = (01, . . . , a n ). 

The following lemma rephrases Lemma [TTI 

Lemma 14: Suppose the quadratic form (1, — a, ja, — 6) is anisotropic over F and (7, — 9) is anisotropic over 
K. Then nonzero elements of A have nonzero determinant. 

Proof: Since (j,—9) is anisotropic over K, there are no nontrivial solutions to the equation 7a; 2 — 9y 2 = 0, 
which implies that 7 7^ 9 (mod A" x2 ), giving us the second condition of Lemma [TTI Next we verify that the first 
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condition of Lemma [TT] is satisfied. We will show that for any u £ K,v <S \/aF, we have 9 ^ ua(u) + 7U 2 . 
Suppose, for the sake of contradiction, that it is not so. Rewriting with u = uq + s/aui, v = y/avi we obtain 
ua(u) + 7i' 2 = Uq — au\ + ^avf — 8 = 0. But this implies that the quadratic form (1, —a, 7a, —6) has a nontrivial 
solution in F, contrary to our assumption. ■ 
The next result is an application of Springer's Theorem. 

Corollary 3: Let Of denote the ring of algebraic integers of F. Suppose a generates a prime ideal of Of and 
that 9,j are non-square modulo aOj?. Then the quadratic form (1,-0,70, —9) is anisotropic. 

Proof: This is an application of Springer's Theorem (Theorem 6.1 lTT9l ). from which it follows that the 
quadratic form (1, — 0,7a, — 9) is anisotropic if and only if the residue forms (1,-6*), (7,-1) are anisotropic 
over the residue field Of/ciOf- Since 9,j are both non-squares (mod clOf), there are no nonzero solutions to 
the equations y 2 = 9z 2 (mod a), y 2 = 7Z 2 (mod clOf)- In other words the residue forms (1,-6*) , (1,-7) are 
anisotropic and we obtain the conclusion. For a more detailed version of a nearly identical proof see Theorem 7.1 

ED. ■ 

Putting the above together we obtain the following. 

Corollary 4: Let a generate a prime ideal of Of- Suppose #,7 are both non-square (mod a) and also suppose 
(7, —8) is anisotropic over K. Then nonzero elements of A have a nonzero determinant. 

Proof: By the previous corollary and Lemma [14] the conditions of Lemma [TT] are satisfied. ■ 

Example 3: Let F = Q(i) and let Q = (3,7)f, so that K = F(y/S). We would like to pick 7, 9 e F which 
satisfy the criteria of the previous corollary, and hence Lemma [TT] Note that 3Z[i] is a prime ideal of Z[i] of norm 
9, i.e., Z[i]/3Z[i] is a finite field Fg of 9 elements. Let 7 = 1 + i, 6 = 1 — i. We can directly verify that both 6, 7 
are non-square modulo 3Z[i] and moreover 7/6* = i ^ F(v / 3) x2 - Hence, the code resulting from ae(Q, Q) is fully 
diverse. 

IV. Examples of MIDO Code Constructions 
In this section we illustrate the general iterated construction using a few well known codes. 

A. An Iterated MIDO Silver code 

The Silver code, discovered in lETI . and rediscovered in ll22l . is given by codewords of the form 



xi —x\ 


+ 


1 




Zl 


* 


x 2 x\ 




-1 




Z-2 


Z *l _ 



where 



Zl 


1 


l + l 


-l + 2i 




%3 


Z2 


~77 


1 + 2i 


1 - i 




X4 



and xi,X2,X3,X4 £ are the information symbols. Consider the number field F = Q(y/—7) and its Galois 
extension K = F(i), with Galois automorphism a : a + bi i-> a — hi, where a,b € F. Note that a is not the 
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complex conjugation since it fixes \/—7. It was shown in |[T4ll that Silver codewords can alternatively be viewed 
as (scaled) matrices over K of the form 

c — cr(d) 
d tr(c) 

where c, d G Z[z] © Z^^ii^-^). To be more specific, let Q = (— 1, — 1)f be a quaternion algebra with F = 
Q(\/— 7). The Silver code is obtained via matrix representation of the elements of the order Ok © jCif of Q, 
that is, for each element c + j<i, with c,d G Q(V — 7, i) we have 

c + ja i-> 

<i ct(c) 

Now for any 9 G F, we have the following mapping, which sends two elements of the Silver code to an element 
of the iterated Silver code: 



ae 



c -a(d) 




e -<r(f) 


d cr(c) 




J o"(e),_ 



(22) 



c -o-(d) 6»CT(e) -6»/ 
d <t(c) 8a(f) 9e 

/ tr(e) o-(d) c 

When 9 satisfies conditions 1), 2) of Lemma [TT] the resulting code will be fully diverse. Next we show this is the 
case for = —17. 

Clearly, 6J—1 = 17 is not a square in K. 

Condition 1) amounts to showing that 6 ^ ua(u) — 71>ct(u), where u G K,v € iF. Writing u = Uo + iui,v = iv\ 
where uq, Ui,v\ G F, and recalling that a : i i-> — i and 7 = —1, we thus want to show that 9 ^ Uq + u\ + v\ 
for any u§,u\,v\ G F. We show that 9 is not a sum of three squares by considering reduction modulo the ideal 
I = ( 1+ f^ ) 3 F - First note that O f /( 1+ ^ =7 ) ^ Z/2Z, O f /I 3 ^ Z/8Z and observe that 9 = -1 (mod 7). 
Using the isomorphism Of/ I 3 — Z/8Z, we verify that —1 (mod /) is not a sum of three squares, and hence 
neither is 9. 

We conclude that the iterated Silver code is fully diverse, in fact it has NVD since F is imaginary quadratic. 

As an example of a basis of an iterated MIDO Silver code, fix the natural basis of the Silver code, and consider 
its image under ae : 

( 1 \ / i \ 

(23) 




l 











•) 


— i 






i 




( 








i 





l + 2i 

-l + i 



.0 
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Otg 0, 




ag 



-2-i 

2-i 

We start with a general lemma which holds for any choice of 9. We use the iterated MIDO Silver code as an 
illustration. However, the same result will hold for iterated codes obtained via ag from Q = (a, —1)f, where a < 
and F is a quadratic extension of Q. In the case of the Silver code F = 7) and a = — 1. 

Lemma 15: The complexity of an iterated MIDO Silver code is at most Od^l 13 ), no matter the choice of 9. 
Proof: We may choose a basis for the Silver code so that it contains the matrices 



1 







i 




-1 




i 





1 


3 


-i 


1 


1 




i 



so that 



Bi = ag 



B 3 = ag 




, , B 2 = ag 



, , B 4 = ao 



i 
-i 

% 

1 



.0 



,0 



are four basis matrices for an iterated MIDO Silver code, independently of 9. It is a direct computation to check 
that 

B k B\ + BiB* k = 0, for all k ^ I G {1, . . . , 4}. 
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M 



(24) 



By choosing an ordering on the basis of an iterated MIDO Silver code such that the first four elements are 
Bi, B2, B3, B4, the above computations show that the associated matrix M, defined in (O, has the structure 

A Mi 
M 2 M 3 

where A is a 4 x 4 diagonal matrix. Hence, by [3] Lemma 2], the decoding complexity reduces from 0(151 16 ) to 

0(\sn. m 

When 6 = Q9' for ( 4 = 1 we get a better reduction in decoding complexity using the map d := a^^g as defined 



in 



•■( 


c -<x(d) 
d <r(c), 


7 


f <r(e), _ 


) 


c 


-a{d) 


c- 






d 


(7(C) 




/&<tU) (V&e 








a(c) 


-d 


VWf 


VWa{e) 




a(d) 


c 



(25) 



We demonstrate this further reduction in complexity in the following lemma. 

Lemma 16: The complexity of the iterated MIDO Silver code (|25j with 6 = -0' is 0(\S\ W ). 

Proof: We show that in fact the code is conditionally 4-group decodable. Using the basis for the Silver Code 
given in (l23l l. in this exact order, define the partitions 

r x = {i, 11}, r 2 = {3, 9},r 3 = {4, io},r 4 = {2, 12}, 

T c = {5,6,7,8,13,14,15,16}. 

It can be verified by direct computation that the iterated MIDO Silver code with 9 = —6' has matrix M, whose 
first 8 rows have the form 

'toootttt 
otootttt 
oototttt 
ooottttt 

where denotes a 2 x 2 matrix and t any 2x2 nonzero matrix. It can thus be seen that the iterated Silver code is 
conditionally 4-group decodable: the four groups of two symbols are clearly distinct. The meaning of conditionally 
group decodable is also readily understood: conditioned on decoding the last 8 information symbols (for a cost of 
Od^l 8 )), the other half of the symbols is decodable by pair (for a cost of OdS*! 2 )) independently, resulting in a 
complexity of 0{ \S\ 10 ). ■ 
Furthermore, for the sake of experiment, we tried the choice of 9 = i, another root of unity which seems very 
natural, even though it belongs to K and thus is not included in the description of our iterated code construction. 
An explicit computation of the M matrix when 6 = i shows that it is of the form ([24-b . Consequently, the decoding 
complexity of the code (F22l when = i is 0(|5| 13 ). 
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M-4,N.2,4-QAM 



M.4,N-2,4-QAM 





10 11 12 



Fig. 4. Comparison among codes with decoding complexity 0(\S\ 12 ) and 0(|S| 13 ) on the left, and of decoding complexity 0(|5| 10 ) on 
the right. 



The performance of two versions of the iterated MIDO silver code, one with 6 = — 1 of complexity OdS| 10 ) and 
the other one with 8 = i and complexity OdS'l 13 ), is illustrated on Fig. [4] On both figures, the x-axis corresponds 
to the signal-to-noise ration (SNR) in dB, and the y-axis to the Block Error Rate (BLER). The signal constellation 
S used is a 4-QAM constellation. On the left of Fig. |U codes with decoding complexity OdS'l 12 ) and 0(|S| 13 ) 
are displayed: we see that the difference among the codes is tiny up to 18 dB. The BHV code [2| has complexity 
0(|S| 12 ) and is known not to be fully diverse. It has a very good performance, and this is only starting from 18 
dB that the loss in diversity starts to be noticed. The SPCOM code ll23l has worse decoding complexity than the 
BHV code, but is fully diverse, and thus gives a way of comparing the gain in diversity with respect to the BHV 
code, which is best seen at higher SNR. The iterated MIDO Silver code with = i also starts to exhibit a loss in 
diversity around 18 dB. The iterated MIDO Silver code with 6 = —1, which is not fully-diverse, is also shown 
for the sake of comparison. The fact that this code is not fully diverse can be seen much earlier than for the BHV 
code, though its final performance in this range of SNR is quite similar. To get a fairer comparison, the iterated 
MIDO Silver code which has decoding complexity 0(|S| 10 ) is compared with two codes of same complexity, based 
on crossed product algebras J5]. The code O5 is not fully diverse, but thanks to a good shaping, has very good 
performance from low to medium SNR. The code O4 is fully diverse, but its performance suffers from the lack of 
cubic shaping. The newly proposed code behaves closely to the code O5. 



B. An Iterated MIDO Golden Code 

Consider the quaternion algebra Q = (5, iWj), with K = Q(i, v5), o : v5 ^ — v5. Let <fi = (1 + v5)/2 denote 
the Golden number, with <j{4>) = (1 — v5)/2. Set j3 = 1 + ia{4>), with <r(/3) = 1 + i(f>. We can view elements of 
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the order {(3q, q 6 Q} of Q as 2 x 2 matrices over via representation in M nxn (K) as 

/3c icr(^d) 
/3d cr(/3c) 

which, up to a scaling by l/\/E, are codewords of the Golden code |[l"3l . ||241 . Now define, for a suitable 9 € Q(i), 
the map 



By Lemma [JJ the determinant of each codeword is in Q(i), moreover it is indeed in since the coefficients 
of each codeword come from the ring of algebraic integers of K. Hence in order to show that the iterated Golden 
code has NVD, it suffices to pick suitable 9, so that the determinant of each nonzero codeword is nonzero. In the 
following lemma we do so by applying Lemma [T4l 

Lemma 17: Let 6 = 1 — i. Then the iterated Golden code arising from ag(Q, Q) is fully-diverse and has NVD. 



Lemma [T4l says that it suffices to pick 9 E Q(i) so that the quadratic form (1,-5,57,-0) i s anisotropic over 
Q(i), and (7,-6*) is anisotropic over K. Recalling that 7 = i, we need to show that (1,— 5,5i,— (1 — i)) is 
anisotropic. Note that 5 factors into 5 = (2 + i)(2 — i) in Q(i), each of the factors generating a prime ideal. By 
Springer's Theorem [19. Theorem 6.1], it suffices to show that the quadratic forms (1, —6) and (—(2 — i), (2 — 
are anisotropic over the finite field Z[i]/(2 + i). Showing that (1, —9) is anisotropic is equivalent to establishing 
that 9 is not a square in Z[i]/(2 + i) = F5. By our choice 9 = 1 — i = 3 mod (2 + i), which is not in Fg 2 = {±1}. 
Similarly showing that (—(2 — i), (2 — is anisotropic over (2 + i) is equivalent to establishing that i is not 
a square in Z[i]/(2 + i), which follows from the fact that 8 does not divide the order of Fg . 

It remains to show that 9/-f $ K x2 . Let z = 9/j = i-p = — i— 1. Suppose z = x 2 . Then zz* = (— i— — 1) = 
2 = (xx*) 2 . However, clearly V2 ^ K, which completes the argument. 

■ 

C. An Iterated MIDO Alamouti Code 

We consider the iterated construction of ExampleQ] with Q = ( — 1, — 1)q, F = Q, K = Q(i) and for x <E Q(i), 
a : x 1— > x* denotes the complex conjugation. We view elements of Q as 2 x 2 matrices over K = Q(i) of the 




/3a ir(/3)r(6) 9t{/3)t{c) i9f3d 

/3b r(/3)r(a) 6»r(/3)r(d) 9/3c 

(3c iT{(3) T {d) r(/3)r(a) i/36 

/3d t((3)t(c) r((3)r(b) (3a 



Proof: 
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form (ETl i. corresponding to the celebrated Alamouti code Q. Now 



ag : 



( 


/ -9* 




c -d* 




9 r _ 




d c* 



H->- 



./• 


-.9* 


0c* 


-ed 


.9 


r 


6»rf* 


9c 


c 


-d* 


/* 


-9 


<7 


c* 


9* 


f 



The case 9 = — 1 corresponds to the quasi-orthogonal code proposed by Jafarkhani 12511 . However, this code is not 
fully diverse, as already noticed in 11251 . This can be remedied by replacing —1 with any negative 9 to satisfy the 
criteria of Lemma [TT] as long as 9/-f = — 9 $ <Q(i) x2 , e.g. 9 = —3. However this code is not full-rate, as it can 
transmit up to 4 complex, equivalently 8 real information symbols. 

In order to increase the rate of such a code, we could instead try to consider the quaternion algebra Q = 
(— 1,— 1)q(^/5), or more generally Q = (— a, — T^yg), with 0,6,7 > 0> as proposed in |26l . This in fact does 
give very good orthogonality relations among the basis elements, as shown in ||26l . However, choosing as a base 
field a totally real quadratic number field means that the basis matrices {D\, 
linearly dependent. To see that, recall that {D\, . . . , Dg} are 



, D s } generating Q over Q will be 





1 







Vb 




\J —a 





D 1 = 






,D 2 = 


,D S = 











1 




Vb 







-a 








—a^/b 





,D 5 = 





-77 













-y/—ay/b~, 


v/7 








A; 



,D 7 



77/ 



-VlVb 
VlVb 

yfjy/^Vb 
FfV^Vb 

We can observe that D\ and D2 are real multiples of each other, implying that 

h 



Si 



a (D u O) 



7 2 



So 



(D 2 ,0) = Vb 



h 
h 



are real multiples of each other as well. This affects the rank of the matrix B in (0, thus that of R; more specifically, 
the matrix R is no longer full-rank. As a consequence, we cannot draw conclusion about fast-decodability of the 
code, despite the sparsity of M. 

This can be avoided by using a totally imaginary quadratic field as base field. The complexity induced in such 
a case is similar to that of the Silver code, discussed above. 



V. Examples of Higher Dimensional Code Constructions 

We now explain how our general construction applies to the case of an asymmetric MIMO channel with 6 
transmit and 3 receive antennas 11271 . 
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C 2 xC 3 



F 3 7 



Fig. 5. Field extension used for the proposed iterative code construction. 



A. Cyclic Algebras of Degree 3 

Let L/F be a Galois extension of number fields with Galois group Ci x C3 ~ Cg, where a denotes the generator 
of the cyclic group C3 and r that of C2 (see Fig. 0. Let A' be the fixed field of a, so that L/K forms a cyclic 
extension of degree 3. Let T> be the corresponding cyclic algebra of degree 3. The representation A of elements 
of T> is given by 

a 7ct(c) ja 2 (b) 
b a(a) 7tx 2 (c) 
c a(b) cr 2 (a) 

where 7 G A is chosen so that X> is division, which guarantees full diversity (|6j. 
Given an element # of K, let ag : T> x 2? — >• Mzft) be the map defined in (TTOb by 

a; 6*r(j/) 

a e : (x,y) 

?7 r(x) 

where x,y are identified with their representation X(x),X(y), and r(x) is defined componentwise, so that 



a, b, c £ L, 



(26) 



a 7<t(c) 7<t 2 (&) 
6 cr(a) 7<7 2 (c) 
c <x(6) c 2 (a) 



a' 7<t(c') 7<t 2 (6') 
6' a(a') 7 cj 2 (c') 
c' cr(6') a 2 (a') 



\ 



(27) 



a 7<t(c) 7cr 2 (fe) 6>r(a') djT<r{c') 6jTcr 2 {b') 
b o-(a) 7ct 2 (c) 0t(6') 0tct(o') 6<7tgt 2 (c') 
c cr(6) cr 2 (a) 0r(c') 0to-(6') 6»rcr 2 (a') 
a' 7<t(c') j<r 2 (b') r(a) 7rcr(c) jTa 2 (b) 
V a(a') 7 ct 2 (c') r(6) ra(a) 7 ra 2 (c) 
c' ct(&') cr 2 (a') r(c) rtr(6) ra 2 (a) 
We now describe one scenario when the image A of ag(T> : T>) satisfies the full diversity property. 

Lemma 18: Suppose that 7 is chosen in F, so that t(j) = 7. Let 6 be an element of K such that t(# 3 ) 7^ 6* 3 
Then A satisfies the full diversity property. 
Proof: 

By Lemma [7] the claim is equivalent to showing that 
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L = Q(( 7 ,i) 



K = Q(i, V7)3 = iV7 



(i) 3 7 = 1 + i 



Fig. 6. Field extension used for the first code example. 



xy- 1 t{xy- 1 ) ^ eh 

Let Z = XY^ 1 and write it as 

u 7<t(«;) 7tJ 2 (z;) 
•2'= y cr(w) 7<j 2 (u>) , u,v,w e L. 
w cr(v) f 2 (u) 
For the sake of contradiction, suppose Zt(Z) = 8I3. Then 

dct(Z)det(r(Z)) = dct(Z)r(dct(Z)) = (9 3 . 

Since det(Z)r(det(Z)) is fixed by r, it should be that 6> 3 is also fixed by t, a contradiction. ■ 
Remark 19: By choosing (9 not in F, A is not multiplicatively closed, thus it does not form an algebra (see 
Lemma [6]). 

We now provide examples of explicit code constructions that fit the above setting. 



B. Examples 

Example 4: T> = (L/K, <T2, 1 + i), with L = Q(C7, i), where £7 is a primitive 7 th root of unity, (J2 '■ Cj ^ (7 
and K is the fixed field L^ 2 ^ . To define the map a$, we consider t to be the generator of Gal (Q(\/—7, i)/Q(i) 
and values f? = \/— 7, f? = — 1. 

The Galois group Gal(L/Q(i)) is cyclic of order 6, generated by the automorphism 

<?3 ■■ (7 C 7 3 - 

Note that the automorphism 

^2 : C7 H' C7 

is of order 3 and hence has a fixed field K which is of degree 2 over Q(i). Therefore K = Q(v^7, i) = Q(V7, i), 
with the Galois (Q(i)-automorphism 

t : V7^-V7, 
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t : 1 1 — y i. 

Claim 1: The algebra T> = (L/K,a2, 1 + i) is division. 
Proof: Note that the ideal 7Z remains inert in Z[i], in other words, its residue field Z[i]/7Z[i] is the finite 
field of 49 elements. We further note that the ideal 7Z[i] is completely ramified in Q((,7,i), hence its residue 
field remains to be of size 49 in K. By local class field theory, the reciprocity map induces an epimorphism 
■^49 — > C3 = Gal(L/K), whose kernel is a subgroup of F^ 9 of order 16, which contains the image of norms of L/K 
(see l20l for more detailed treatment of this type of argument). Hence if an element is a norm of Q(C7, i)/Q(i, V7), 
then its image in F|g is an element whose order divides 16. 

This allows us to conclude that 7 = 1 + i, which has order 24 in F£ 9 is not a norm in the extension L/K, hence 
the cyclic algebra (L/A', o^, 7) is division. This completes the claim. ■ 
The choice 7 = 1 + i satisfies the condition t(j) = 7, as required by the iterated construction. Note that 
6 = iO' = satisfies the condition t(0 3 ) ^ 6 3 of Lemma [T8l Since = Q(i), Lemma [7] implies that the 
resulting code has NVD. From Lemma [TO] the iterated code obtained via a iv ^(X>, £>) defined in (fT4b will hence 
inherit fast-decodable properties of T>. Next we fix a convenient basis for the code arising from T>. This code 
should be of Q-rank 18, so that the iterated code is of desired Q-rank 36. Let {v\, z^, ^3} be a A"-basis of L, such 
e.g. {1, ( 7 + £7 1 , Q + (7 2 }. The representation A of V gives a family of 3 x 3 matrices of the 



that all Vi £ 
form 



2 

E 

k=0 



c k 
tr(c fc ) 









where 



7 

1 
1 

and Co, c%, C2 are A"-linear combinations of {v\, u^, ^3}, that is 

2 3 - 



fc=0 j=l 



a{vj) 
a\v 3 )_ 

We can thus encode 9 elements of K. Since K is of degree 2 over Q(i), this is more rate than we can support, 
namely 18 complex symbols instead of 9. Let us write d £ K as 

d = do + d\i + yJdi + v7id3, di,d2,d%,d± £ Q. 

By setting two dt to zero, we will obtain the right rate. Instead of the Q-basis \/7,i\/7} of K, we will thus 
use a subset {/^i,/^}- For reasons of fast-decodability we insist that be real (pi is 1 or vT) and fi2 be totally 
imaginary (^2 is i or i\7). We obtain from a A-basis of T> the following Q-linearly independent set for j = 1,2,3: 



Dj = (ixVj, D 3+j = fx 2 Vj. 
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D 6+j = niVjT, D 9+j = fi 2 VjT, 



D 12+j = /iiKT 2 , D 



15+j 



where we denote 









trfa) 
a 2 { Vj ) 



for short. We have 



HiVj^ + ^VkvlV* = M + to)VjVk = o 

using the fact that Vj, Vk have only real coefficients, and commute since they are diagonal. Furthermore, \i\ was 
also real, while ji 2 is totally imaginary. Next 







(n(j4 ^ 


-^WjTrVk = 














7 




1 




l7| 2 








rr* = 


i 







1 







1 










1 




7* 










1 



The same computation holds for the last set of basis elements, involving the generator T 2 . 

Mi^r^r 2 *^* + fi 2 v k r 2 ^y 2 *v* = o 

At the price of sacrificing full-diversity, we may choose 9 = — 1. This will result in additional relations: 

a e (D i ,0)a e (0,D j y + a e {0,D j )a e (D j ,0)* = 

for all i, j = 1, . . . , 6, implying decoding complexity Od^l 30 ) using |3] Lemma 2], rather than the worst case 
0(\S\ 36 ) for the given code. 

Example 5: V = {L/K, a%, 3) with L = (QKC7), cr 2 : C7 ^ C? an d K = Q(V~ 7)- To define the map a$, we use 
t : y— 7 M> — \/— 7 and consider 6* = \/— 7. 

Claim 2: The algebra 2? = (L/K,a2,3) is division. 
To see this, we use similar arguments as in Example |4] to demonstrate that 3 is a non-norm element in the extension 
L /K. Namely, 3 multiplicatively generates nonzero elements of F7, which is the residue field of the ideal above 7 
in K. Therefore the image of 3 in F^ cannot lie in the kernel of the epimorphism Fj — s- C3. Hence the algebra 
T> = (L/K, (72, 3) is division. The representation A of T> gives a family of 3 x 3 matrices of the form 





EE** 

k=0j=l 



a(uj) 
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L = 




K = Q(iV7) 3 = iV7 



F = Q 9 7 = 3 



Fig. 7. Field extension used for the second code example. 



where cjy £ K and {ui, V2, ^3} is a A'-basis of L. Each matrix encodes 9 elements of A, or 18 elements of Q, 
since this time A is of degree 2 over Q. For the iterative construction, first we observe that 6 = i\J~7 is an element 
whose cube is not fixed by r. Hence, by Lemma Q~8] the elements in the image of <x iy m have nonzero determinant, 
moreover since = Q, Lemma [7] implies that the resulting code has NVD. We set the code to be the matrices 

in the image of <S iv /g(2?, 2?), and by Lemma [TO] fast-decodable properties of T> are inherited. Similar computations 
as in the previous example can be done. 

VI. Conclusion and Future Work 

Motivated by the code design of fast decodable space-time codes for asymmetric MIMO channels, we proposed 
an iterated code construction, that maps a space-time block code coming from a cyclic division algebra into a 
new space-time code of twice the original dimension. We studied the algebraic structure of the newly obtained 
codes, which are actually forming interesting families of division algebras, whose center and maximal subfield are 
computed. In particular, we gave a condition for full diversity. Particular attention was paid to the quaternion case, 
corresponding to the MIDO case. 

It is surely possible to obtain other (possibly better) codes from the general construction that we proposed. It 
would be a valuable contribution to establish a bound on the best decoding complexity that codes coming from 
division algebras can have. From a purely mathematical point of view, it would be of interest classify the division 
algebras that are built through the proposed iterative construction. 
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